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ABSTRACT 

The initial conditions of massive star and star cluster formation are expected 
to be cold, dense and high column density regions of the interstellar medium, 
which can reveal themselves via near, mid and even far-infrared absorption as 
Infrared Dark Clouds (IRDCs). Elucidating the dynamical state of IRDCs thus 
constrains theoretical models of these complex processes. In particular, it is 
important to assess whether IRDCs have reached virial equilibrium, where the 
internal pressure balances that due to the self-gravitating weight of the cloud 
plus the pressure of the external environmental. We study this question for the 
filamentary IRDC G035. 39-00. 33 by deriving mass from combined NIR & MIR 
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extinction maps and velocity dispersion from C 18 (1-0) & (2-1) line emission. In 
contrast to our previous moderately super-virial results based on 13 CO emission 
and MIR-only extinction mapping, with improved mass measurements we now 
find that the filament is consistent with being in virial equilibrium, at least in 
its central parsec-wide region where ~ 1000 M© snakes along several parsecs. 
This equilbrium state does not require large-scale net support or confinement by 
magnetic fields. 

Subject headings: ISM: clouds, dust, extinction — stars: formation 

1. Introduction 

Identified by obscuration of the mid-infrared (MIR) (i.e. ~ 10/im) Galactic background, 
Infrared Dark Clouds (IRDCs) are likely to be representative of the initial conditions of mas- 
sive star and star cluster formation, since their high mass surface densities (E > O.lgcm -2 ) 
and densities («h > 10 4 cm 3 ) are similar to regions with such star formation activity (e.g. 
Teyssier et al. 2002; Rathborne et al. 2006; Tan 2007; Butler & Tan 2009, hereafter BT09; 
Zhang et al. 2009; Ragan et al. 2009; Butler & Tan 2012, hereafter BT12). 

The kinematics of IRDCs can be measured via their molecular line emission to determine 
if they are gravitationally bound and/or in virial equilibrium. Hernandez & Tan (2011, 
hereafter HT11) used 13 CO(1-0) emission from the Galactic Ring Survey (Jackson et al. 
2006) to measure velocity dispersions in two filamentary IRDCs (F & H in the BT09 sample). 
£ was estimated by averaging two methods: (1) MIR extinction (MIREX) mapping (BT09), 
given an assumed MIR opacity per gas mass and a foreground correction based on an analytic 
model of Galactic hot dust emission; (2) 13 CO emission, given an assumed abundance of this 
isotopologue. Cloud mass was then calculated assuming near kinematic distances to the 
IRDCs. A filamentary virial analysis following Fiege & Pudritz (2000, hereafter FP00) was 
performed in several orthogonal strips along the filaments. HT11 concluded surface pressure 
terms are dynamically important, suggesting that the filaments have not yet reached virial 
equilibrium. 

Here we revisit this analysis for one of the filaments, G035.30-00.33 (H in the BT09/BT12 
sample), which is 2.9 kpc distant with a few thousand solar masses of material spread over 
about 4 pc (projected) at its northern end. Although there are several 24 /iin sources seen 
towards the filament (Carey et al. 2009), which are likely to be embedded protostars, most 
of the region appears MIR dark and starless. For our kinematic measurements, we use 
IRAM 30m observations of C 18 (1-0) and (2-1). The results of these and other molecular 
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line observations are being presented in a series of papers on the formation and evolution 
of this filamentary IRDC Paper I (Jimenez-Serra et al. 2010) presented maps of SiO, CO, 
13 CO and C 18 0. Widespread SiO emission was observed, perhaps suggesting the presence of 
large-scale shocks that may have been involved in forming the filament. Paper II (Hernandez 
et al. 2011) compared the C 18 and BT12 MIREX maps, showing CO suffers widespread 
gas phase depletion in the IRDC by factors of up to f' D ~ 5, i.e. just 1 in 5 CO molecules 
remain in the gas phase. Paper III (this Letter) performs a filamentary virial analysis of the 
IRDC. Paper IV (Henshaw et al. 2012) studies the detailed kinematics of the filament and 
its surroundings, including analysis of the dense gas tracer N2H + . 



2. Mass Surface Density from Extinction Mapping 

We use the MIREX E maps of BT12 combined with near infrared (NIR) (i.e. J,H,K 
[UKIRT Infrared Deep Sky Survey, Lawrence et al. 2007]) extinction maps (Kainulainen et 
al. 2011a) to yield a final map that has been presented by Kainulainen & Tan (2012, hereafter 
KT12). The MIREX map is derived from the 2" resolution Spitzer IRAC 8 /mi (Galactic 
Legacy Mid-Plane Survey Extraordinaire [GLIMPSE]; Benjamin et al. 2003) image (Figs. 
UK & b). A MIR (band and background spectrum weighted) dust opacity per unit total 
mass of fi^m = 7.5 cm 2 g^ 1 was adopted, consistent with the moderately-coagulated thin ice 
mantle dust model of Ossenkopf & Henning (1994, hereafter OH94) and a gas-to- (refractory- 
component-) dust mass ratio of 156 (Draine & Lee 1984). We estimate that the uncertainties 
associated with these choices are ~30%, including those due to the grain composition and 
size distribution (e.g. the OH94 uncoagulated thin ice mantle model has opacities 0.83 times 
smaller; the maximally coagulated model 1.11 times larger; the bare grain, uncoagulated 
model 0.67 times smaller) and the gas-to-dust mass ratio (e.g. Draine (2011) estimates a 
ratio of 141 from depletion studies, 0.90 times smaller than our adopted value). Bright MIR 
sources lead to "holes" in the MIREX map. These cover just a small fraction of the projected 
area of the IRDC. They are excluded from the calculation of mean E. 

The effect of BT12's use of an empirical foreground estimate (from observed saturation 
in independent cores) compared to BT09's analytic estimate is to increase E, especially in 
the densest regions, where it rises by up to a factor of ~3. The average increase of E inside 
the elliptical area of the IRDC defined by Simon et al. (2006) is a factor of 1.8. 

In regions of low E near the edge of the cloud, the MIREX maps underestimate E since 
the technique is based on relative extinction of dense regions compared to the surroundings. 
This is now corrected for by combining the MIREX maps with NIR extinction maps (KT12). 
In regions where the NIR extinction map has Ay < 10 mag (~ 0.04 g cm -2 ) it is used to set 
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Fig. 1— IRDC H (G035.30-00.33). Left to Right — (a): Spitzer GLIMPSE IRAC 8 image. 
Following HT11, four E-W strips (1-4 from N to S) have been drawn across the filament. Red 
dashed vertical lines delineate either side of the "Outer Filament". Magenta solid vertical lines 
define the "Inner Filament" regions, each half the width of the Outer Filament and centered on the 
center of mass of the strip. The image has 1.2" pixels and PSF with 2" FWHM. (b): Total mass 
surface density, E, derived from (a) using the MIREX mapping method (BT12) with low mass 
surface density normalization set by NIR extinction mapping (KT12). (c): C 18 0(2-l) emission 
integrated over velocities from 40 — 50 km s^ 1 (HT11), a pixel scale of 5", and a beam size of 11". 
(d): Total ID velocity dispersion, a, derived from the C 18 0(2-l) 40 — 50km s _1 spectra, (e): 
Normalized CO depletion factor, f' D , for regions with Scrso > 0.005 gem -2 (see text). 

the zero point of the MIREX map. In higher column density regions, only the MIREX map 
is used, now with the zero point interpolated from the surroundings. 
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3. C 18 Observations to estimate Velocity Dispersion and CO Depletion 

The C 18 J = 1 — and 2 — 1 transitions were mapped (On-The-Fly) over a 2' x 4' 
area with central position a(J2000) = 18 h 57 m 08 s , 5(J2000) = 02°10 / 30" (I = 35.517°, 
b = —0.274°) using the IRAM (Insituto de Radioastronmfa Millimetrica) 30m telescope 
in Pico Veleta, Spain in August & December 2008 with a combination of the ABCD (for 
C 18 O(l-0)) and HERA (C 18 0(2-l)) receivers (see Paper II for more details). At ~ 110 GHz 
(J = 1 _> 0), the beam size is 22", while at ~ 220 GHz (J = 2 -> 1) it is 11". Data 
were calibrated using GILDAS (CLIC). Spectra were regridded to a velocity resolution of 
0.2 km s _1 . 

The C 18 emission from the IRDC is between velocities 40 — 50km s^ 1 (HT11). The 
C 18 0(2-l) integrated intensity map is shown in Figure [lb. Standard analysis methods (cor- 
recting for excitation temperature and optical depth variations across the position-velocity 
cube) were used to derive the distribution of CO-emitting gas as a function of velocity (Pa- 
per II) and thus a map of total (thermal + nonthermal) ID velocity dispersion, a, assuming 
mean particle mass /i = 2.33m p (i.e. %e = 0.2nn 2 ) and temperature of 15 K (Figure [Tbi). 
We assume a is measured to about 10% accuracy. Note that the IRDC is clearly revealed as 
a region of relatively low a. Mass surface density assuming a constant C 18 abundance was 
derived via S C180 = 7.65 x 10~ 2 (iV cl8 o/10 16 cm~ 2 ) g cm -2 (Paper II). 

The ratio of mass surface density from NIR and MIR extinction, E (smoothed to the 
pixel scale of the CO data), and that from assuming a constant C 18 abundance, Eciso 
defines the CO depletion factor f D = E/E C180 (averaged along a sight line through the 
cloud). To examine relative changes in depletion, following Paper II, we define a normalized 
depletion factor f' D = Bf D , where B is a scaling factor so that f' D is on average equal to 1 
in regions with 0.01 < E/g cm~ 2 < 0.03. We find B = 0.71 compared to the Paper II value 
of 1.47 (Case 1 [i.e. no envelope subtraction] HiRes) due to changes in the E map, discussed 
above. For example, if CO is undepleted in this region and E correctly estimated, this would 
imply we underestimated the true C 18 abundance by ~30%. The map of f' D is shown in 
Figure le, showing values for wherever Eciso > 0.005 gem -2 . Given the improved E map 
due to combination with NIR extinction data, this f' D map supercedes that of Paper II: in 
particular it extends to regions of lower mass surface density. However, the basic structure 
is unchanged, with maximum values of f' D ~ 5 in the highest column density cores. 



6 



4. Filamentary Virial Analysis 

4.1. Filament and Envelope Geometry 

We assume the filament is locally cylindrical, with radius Rf, and symmetrically em- 
bedded in a cylindrical envelope of outer radius R e (and inner radius Rf). To derive the 
mean mass surface density of the filament £/, we need to subtract the contribution from the 
overlapping portion of the envelope, S e! on- For a uniform density envelope and for negligible 
net contribution of material outside R e , we first measure £ ei0 ff , i.e. of the observed envelope 
region that is offset in projection from the filament, by averaging the two parts of each strip 
that are at projected distances between Rf and R e from the filament's central axis. From 
simple geometry, the mean mass surface density of the whole envelope, E e , is related to E ei0 ff 
via 

- _ 7rsin 2 fl(l - cosg) - 

he ~ 29-sm29 ^ {L) 
where cos^ = Rf/R e . Thus, since S e = (Rf / R e )T> etOI1 + ([R e — i?j]/i? e )S eoff , we have 

As shown in Figure 1, we consider two cases: (1) An "Outer Filament" (red lines in 
Figured]) of diameter 2.05' in R.A. and for which R e = l.bRf (set by the extent of the C 18 
observations), i.e. 9 = 0.841 radians, S e = 0.845£ eiO g- and S 6i0n = 0.768£ ejO ff. Each of the 
four strips is centered at a(J2000) = 18 /l 57 m 08.02 s and is 1.12' wide in Dec, with the strip 
2-3 boundary at 5( J2000) = 02°10'35.7". This was the same size and location of the filament 
region adopted by HT11. (2) An "Inner Filament" (magenta lines in Figured]) that we set to 
have half the Outer Filament width and to have R e = 2.0Rf, i.e. 9 = 1.047, S e = 0.959E e , o ff 
and E ei on = 0.918S e)O ff . Now within each strip, the location of the Inner Filament is centered 
on the center of mass of the whole strip, so this case follows the structure of the filament 
more accurately. The correction factors for the above two cases were applied for estimating 
S e and E/. This division of the cloud into filament and envelope regions is quite idealized, 
but is necessary for comparison with analytic models. 

The velocity dispersion of the envelope, a e , is measured from the average spectra of the 
E and W envelope "off" regions in each strip. Then the above factors for estimating contri- 
butions from an overlapping envelope are also applied to the C 18 spectra: the spectrum of 
the filament is derived from the observed "on" spectrum minus the envelope "off" spectrum 
scaled by 0.768 and 0.918 for Outer and Inner cases, respectively. This subtracted spectrum 
is used to derive the total ID velocity dispersions in the filament, <jf. For both o e and <jf, 
only those parts of the spectra with signal > 1 standard deviation of the noise level are used. 



7 



For each strip and for both Outer and Inner cases, the properties of the filament and 
envelope are listed in Table [TJ In addition to the 30% error we assume for the estimates of 
the directly observed values of E (see §2}, we also include an additional 30% error (added in 
quadrature) for parameters that depend on the assumed cylindrical filament geometry, i.e. 
S e and £/. Without direct observational constraints, we further assume the filament axis is 
inclined by % = 60° ±15° to our line of sight (90° would be in the plane of the sky). This 
uncertainty contributes up to 22% uncertainty to quantities such as the mass per unit length 
of the filament. The distance to the filament is assumed to be 2.9 ± 0.5 kpc (Simon et al. 
2006). This uncertainty is equivalent to non-circular motions of ~ 8km s _1 . 

4.2. Comparison to Fiege-Pudritz Models of Virial Equilibrium 

We perform a filamentary virial analysis following FP00, who showed pressure-confined, 
non-rotating, self-gravitating, filamentary (i.e. lengths ^> widths) magnetized clouds that 
are in virial equilibrium satisfy 



Here: P e = p e o\ (with p e = /ih^h>, with mass per H nucleus of /in = 1.4m p = 2.34 x 10~ 24 g) 
is the external pressure at the surface of the filament assumed to be equal to the mean 
envelope pressure; Pf = Pf^f is the mean total pressure in the filament; m/ is the mass 
per unit length of the filament; m vir j = 2aj/G is its virial mass per unit length; M.f is its 
magnetic energy per unit length (associated with large-scale static magnetic fields providing 
either support or confinement); and Wf = —m 2 jG is its gravitational energy per unit length. 

We can measure P e /Pf and m//m v ; r j for the Outer and Inner Filament cases for each 
strip (Tabled]). No information is currently available for Aif- In Figure [2] we thus compare 
the observed filament properties to a series of models that vary the contribution of magnetic 
support (-M//|W/| > from poloidally-dominated B- fields) or confinement (A"f//|W/| < 
from toroidally-dominated B-fields). 

Figure [2] indicates that as one zooms from the Outer to Inner Filaments, P e /Pf decreases 
and mj/m vir j tends to increase. Within the uncertainties, the Inner Filament (and much 
of the Outer) is consistent with the model of virial equilibrium not requiring any magnetic 
support or confinement (7W / / 1 1 = 0). Stated another way, for the Inner Filament the 
internal pressures are factors of a few greater than external. Such pressure enhancements 
are consistent with those expected from self-gravity given the observed masses and velocity 
dispersions. Note that the relative positions of the strips in Figure [2] are more accurately 




(3) 
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Table 1. IRDC Filament Properties 



Cloud property > Outer Filament (strips 1 — 4, Total) Inner Filament (strips 1—4, Total) % Error 

(R f = 0.865 pc, R e = 1.5R f ) (R f = 0.432 pc, R e = 2.0R f ) 
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30% 
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2.62 


2.93 


4.25 


3.74 


4.75 


4.60 


4.34 


30% 


E/ + S e ,on (10- 2 gcm- 2 ) 


5.58 


8.12 


7.18 


7.77 


7.17 


7.31 


12.5 


10.3 


11.6 


10.4 


30% 


t e (lO^gcm" 2 ) 


3.09 


2.26 


2.33 


3.03 


2.67 


3.67 


3.32 


3.62 


3.43 


3.52 


42% 


Y. f (lO^gcm" 2 ) 


2.77 


6.01 


5.07 


5.02 


4.74 


3.97 


9.48 


6.98 


8.51 


7.24 


42% 


M e (M ) 


361 


264 
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354 


1250 


287 


259 


283 


268 


1100 


55% 


M f (M ) 


216 
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395 


392 
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272 
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1130 


55% 


n H , e (10 3 cm- 3 ) 


3.27 


2.39 


2.46 


3.21 


2.83 


4.32 


3.91 


4.27 


4.04 


4.14 


51% 


nuj (10 3 cm~ 3 ) 


2.44 


5.35 


4.47 


4.43 


4.18 


7.01 


16.7 


12.3 


15.0 


12.8 


51% 


m f (Mq pc" 1 ) 
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142 


340 


250 
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259 


51% 
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46.14 
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45.80 


45.19 
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20% 
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55% 


Pe/Pf 
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0.499 


41% 


t sJ = 2R f /a f (Myr) 
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known than the error in any one position, since systematic uncertainties are shared. These 
relative positions, especially that of strip 4 with respect to 2 and 3, further support the 
hypothesis that much of the filament is in virial equilibrium with a common, small degree 
of support from large-scale B- fields. 

We have tested the sensitivity of our results to our adopted method of envelope subtrac- 
tion: if no correction to filament properties due to the presence of an overlapping envelope 
is allowed for, then for the Inner Filament we derive <jf to be 12% lower, P e / Pf to be 
23% lower, and mf/m vir j to be 30% higher. These values are still consistent with virial 
equilibrium without needing large-scale magnetic support or confinement. 

Our results are different from those of HT11, who concluded that, for the Outer Fil- 
ament, the surface pressures were all greater than the internal ones. For Outer Fila- 
ment strips 1 to 4, they found P e /Pf = 3.06,1.74,3.32,5.25. Our values are factors of 
0.46, 0.30, 0.21, 0.25 smaller, respectively. For these strips HT11 found m//m vir j = 
0.0789,0.161,0.269,0.367. Our values are factors of 2.1, 2.7, 1.5, 1.4 larger, respectively. 
These differences are mostly due to our use of improved £ estimation methods. The BT12 
MIREX map improves the estimate of MIR foreground intensity, which increases the derived 
values of £ in the denser parts of the filament. Furthermore, we now no longer use CO line 
intensities to derive S: HT11 used 13 CO line intensities as one method and then averaged 
with the BT09 MIREX value. The 13 CO-derived values were generally lower than the BT09 
MIREX values, which we now attribute as due to CO depletion (Paper II). Thus our esti- 
mates of £ and mass have now increased by factors of a few compared to HT11, causing 
P e /Pf to decrease and m//m vir j to increase. Our present study uses C 18 to measure veloc- 
ity dispersions, which is less prone to being affected by optical depth effects than the 13 CO 
used by HT11 (although HT11 did correct for optical depth). For the Outer Filament strips 
1 to 4, HT11 found a f = 1.52, 1.34, 1.03, 0.995 km s _1 . Our C 18 0-derived values are factors 
1.05, 1.10, 1.33, 1.26 times larger. The higher angular resolution of the C 18 0(2-l) observa- 
tions, allow us to probe to smaller scales, i.e. the narrower Inner Filament case, where line 
widths are seen to decrease by factors of about 0.8 compared to the Outer Filament (Figure 
Id). Therefore the overall velocity dispersions adopted here for the filament are quite similar 
to those of HT11 and the differences in our results are thus rather caused by the higher E's 
now derived. 



5. Discussion 



There is some evidence that the IRDC H filament may have formed recently from 
converging flows of molecular gas. Jmienez-Serra et al. (2010, Paper I) reported widespread, 
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more than parsec-scale SiO emission from the filament, which may have resulted from large- 
scale shocks with speeds of at least several km/s. Henshaw et al. (2012, Paper IV) studied 
the kinematics of the region and found evidence for the formation of the main filament 
via the merging flow of surrounding filaments observed in C 18 (with typical densities of 
^H.flow ~ 10 3 cm~ 3 ). They find that the line of sight relative velocity between each component 
of the merging flow and the main filament is t>fl ow ~ 3 km s -1 . 

If the main filament has formed relatively recently, is this consistent with its observed 
state of near virial equilibrium? It should take at least a signal crossing time, t s j = 2Rf/<jf 
for a region of the filament to settle into virial equilibrium. These times are ~ 0.8 Myr for 
the Inner Filament (Tabled]). To form the filament from two converging flows takes a time 
tformj = 0.236(771^100^^(3 pc- 1 )/[(i? / /pc)(i; flow /3 km s- 1 )(n H ,flow/10 3 cm" 3 )] Myr. Applying 
this to the average properties of the Inner Filament (Rf = 0.432 pc, m/ = 259 M Q ), yields 
tformj = 1.4 Myr (see also Paper IV). Thus, even in the scenario of a recently formed filament, 
enough time should have elapsed for it to have settled into virial equilibrium. 

The widespread CO depletion observed in Paper II and here (Figure le), also constrains 
the age of the filament. The CO depletion time is tpj = 0.16(nH i //10 4 cm -3 ) -1 Myr, assum- 
ing a sticking probability of order unity (Tielens & Allamandola 1987). Table [1] lists toj 
for the different strips. For the Inner Filament, these are relatively short, ~ 0.1 Myr, which 
thus provides a lower limit for its age. 

The implications of virial equilibrium of self-gravitating IRDC filaments are profound. 
If this result for G035. 39-00. 33 applies more generally to IRDCs, then it indicates that the 
initial conditions for star formation, including the cores that form massive stars and the 
clumps that form star clusters, are created from environments where approximate pressure 
equilibrium has been established. The value of this pressure is set by the self-gravitating 
weight of the larger scale cloud, i.e. the IRDC, which dominates over the pressure of its 
surrounding environment (unlike for most clumps in GMCs: e.g., Bertoldi & McKee 1992; 
Kainulainen et al. 2011b). This would confirm a basic assumption of the initial conditions 
of the Core Accretion model of massive star formation of McKee & Tan (2002, 2003) and is 
also expected under the scenario of Equilibrium Star Cluster Formation (Tan, Krumholz & 
McKee 2006). The fact that our results differ from those of HT11, highlights the importance 
of improved estimates of masses and surface pressures of IRDCs. These effects may help 
explain other reported discrepancies between dynamical (i.e. virial) masses and true masses 
(e.g. Battersby et al. 2010). Measurement of large scale magnetic field strengths and 
geometries in IRDCs would help to further constrain the models. 

We thank Bruce Draine and an anonymous referee for helpful comments and discussions 
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Fig. 2. — Ratio of surface to internal pressure, P e /Pf, versus ratio of mass per unit length 
to virial mass per unit length, m//m vir /, for strips 1 to 4 of IRDC H. Red line and points 
labelled # Q show the "Outer Filament" . Magenta line and points labelled 7^ show the "Inner 
Filament". The error bars in the lower-left show typical uncertainties. The smooth curves 
show the conditions satisfied by equation ([3]) for confining magnetic fields with A^j/jW 7 /! < 
(dotted lines), supporting magnetic fields with A4f/\Wf\ > (dashed lines), and no magnetic 
fields, i.e. .M//IW 7 /! = (long-dashed line). The results for the Inner Filament indicate a 
dynamical state consistent with magnetically-neutral virial equilibrium. 

that improved the paper. JCT acknowledges support from NSF CAREER grant AST- 
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